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Abstract 

Qh ■ Using the translation method of Tartar, Murat, Lurie, and Cherkaev bounds are derived 

on the volume occupied by an inclusion in a three-dimensional conducting body. They 
assume electrical impedance tomography measurements have been made for three sets of 
pairs of current flux and voltage measurements around the boundary. Additionally the 
conductivity of the inclusion and surrounding medium are assumed to be known. If the 
boundary data (Dirichlet or Neumann) is special, i.e. such that the fields inside the body 
would be uniform were the body homogeneous, then the bounds reduce to those of Milton 
and thus when the volume fraction is small to those of Capdeboscq and Vogelius. 

Keywords: inverse problems, size estimation, electrical impedance tomography 

^ ■ 1 Introduction 

O 

One of the most basic problems in electrical impedance tomography is to estimate the volume 
fraction of an inclusion in a conducting body from one or more sets of electrical measurements 
of current flux and voltage pairs around the boundary. If the boundary data (Dirichlet or 
Neumann) is special, i.e. such that the fields inside the body would be uniform were the body 
homogeneous, then sharp bounds on the volume fraction were derived by Nemat-Nasser and 
Hori [2TJ |U] (they did not present their universal bounds as bounds on the volume fraction, but 
such bounds can be easily derived from their results) and by Milton |18j . These reduce to the 
asymptotic bounds of Capdeboscq and Vogelius [51 [6] in the important case when the volume 
fraction is small. They are the natural generalization of the well-known Hashin-Shtrikman- 
Tartar-Murat-Lurie-Cherkaev bounds [71 l8l l20l 1241 1141 [15] for composites extended to two-phase 
bodies of arbitrary shape. 

For other boundary conditions bounds, using the information from one measurement pair, 
were derived by Kang, Seo and Sheen [13], Ikehata [ID] , and Alessandrini, Morassi, Rosset, 
and Seo [21 El [I]. These other bounds involve constants which are not easy to determine. 
Another approach was taken by Kang, Kim and Milton [H] who recognized that the translation 
method approach of Tartar, Murat, Lurie, and Cherkaev [231 12Q1 12H El E3 EE] could be directly 
applied to bound the volume fraction using several measurement pairs, with arbitrary boundary 
conditions. The constants in the resulting bounds can be easily evaluated from the boundary 
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data. Furthermore the bounds were sharp if the boundary conditions and inclusion were such 
that the fields were uniform in the inclusion, and shapes of inclusions having this property were 
found. For special boundary conditions they recovered the bounds of [IB] and also provided 
new bounds, even for this case. Those new bounds also reduced to the bounds of Capdeboscq 
and Vogelius [SJ E] for asymptotically small volume fractions. The same approach was applied 
to bound the size of an inclusion in a body for two-dimensional elasticity [19] (where another 
method, the method of splitting, was introduced) and for the shallow shell equations [12J. 

In this paper we extend the translation method approach to three-dimensional, two-phase, 
conducting bodies. This is relatively straightforward for one of the bounds, as there is a null- 
Lagrangian associated with the electric fields, but requires additional care in the case of the 
other bound as there is no null-Lagrangian associated with three-dimensional current fields and 
instead one uses quasiconvex functions [23] . 



2 Upper bounds 

2.1 Derivation of the upper bounds 

Consider three potentials V\, V2, V3 which satisfy V • aWi = in O, i = 1, 2, 3, and Vi = V®, 
and define = — VVJ and jj = ae. L . We rewrite this as 

V • crVV = in tt, V = V° on dtt, (2.1) 

where V = (Vj., V2, V3) is a vector potential, and following Tartar [23] we define matrix valued 
fields 

E = -VV, J = o-E, (2.2) 

where the notation is followed that E has ei,e2, e3 as columns (not rows) and J has ji , J2 : J3 as 
columns. 

We desire to use information about q± = — ji • n, = —J2 • n, qs = — j 3 ■ n and V%, V2, V3 
on d£l to generate an upper bound on f\ (the volume fraction of the phase with the highest 
conductivity). The upper bound is expressed in terms of the response (or measurement) matrix 
A defined by 

A := (E T cjE), (2.3) 
which is computable from the boundary data: 

^kb-^^kL-^- (24) 

We suppose that 

(E) = I. (2.5) 

We emphasize that this condition is achieved by taking a linear combination of the boundary 
data since 

<E> = ML- vTn ' <2 ' 6) 

We further assume that the response matrix A is diagonal, namely, 

A = diag[Ai,A 2 ,A 3 ], (2.7) 

by rotating the body and taking a further linear combination of V°. In fact, since A is 
symmetric, there is an orthogonal matrix R such that R r AR is diagonal. We then take 
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R T V°(Rx) as the new Dirichlet data. We emphasize that the condition (|2.5p still holds, since 
(E) transforms to R T (E)R = I. 

We introduce the 4th order tensor T = (^Tijki) such that 

TP = Tr(P)I - P T , (2.8) 

where P is a 3 x 3 matrix, i.e., 

Tijki = Sij5ki — SuSjk- (2.9) 

The tensor T is associated with a null-Lagrangian and has the property that T(VV) is diver- 
gence free, V • T(VV) = 0. In fact, since 

T( VV, i3 - = % g-||, (2.0) 

we have 2 2 

dx^ ^ dxjdxk dxjdxi ^ 

Note that the Einstein convention for summation is being used. If VW and V( 2 ) are two 
solution vectors, then the matrix 

/ (VVW) r T(VV( 2 ') = -L f (V«f nT(VV< 2 )) (2.12) 
l"l Jn p'l Jan 

can be computed from the boundary data as nT(VV^ 2 ^) only depends on the tangential deriva- 
tives of V^ 2 ) on dtl (V^ 2 ) and n are row vectors). To see this it is enough to notice that the 
i-th component of the row vector nT(VV) is 



We also note that 



/ (VV«) T T(VV( 2 )) = -L [ (VV( 2 fT(VVW), (2.14) 
P'l Jn |"| Jn 

since T is symmetric in the sense that Tijki = Ty-uj for all i,j, k, I. 

Consider the quadratic form Tr(K" r E r T(EK)) on 3 x 3 matrices K. There is a unique 
symmetric 4th order tensor, say M = (Mijki), such that 

Tr(K T E T T(EK)) = K : MK, (2.15) 

where A : B = AijBij. M being symmetric means M^ki = M}.uj for all i,j, k, I. Note that 
M ijH is given by 

M^ = Tr<(EI^) T T(EI fc ')), (2.16) 

where V J is the elementary matrix whose entry is 1 and other entries are 0. We emphasize 
that Mijt-i can be computed by the boundary data as was seen in (|2,12p . In fact, if we let V 1 - 7 be 
the vector such that its jth component is Vi and other components are 0, then (VV)! 1 - 7 = VV 9 . 
Thus we have from (pTH|) and (j!TT3l) that 



I ((VV)F) T T((VV)I H ) = -L f (VV^') r T(VV H ) 
" Jn " Jn 



1 f dVk {Y ij )Tn _ (yij )T n( y yfcZ ST _ 

l"l Jen 5x z 
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Therefore we get 



ijkl 



1 



Vi 



an 



dVk dV k 
dxi 



ni 



We may rewrite this in a symmetric form: 



Mijki = 



1 

W\ 



diii duk dui duk 
dxj dxi dx\ dx 



j J 



(2.17) 



(2.18) 



where ui is the solution to Aui = in Q and ui = V® on dQ, i = 1,2,3. From this one can 
immediately see that M is symmetric, 



if % = k or j = I, 



(2.19) 



and 



M ijk i = -Mukj- (2.20) 

If affine boundary data are used, namely, V® = Xi for i = 1,2, 3, then one can easily see that 
M1122 = M1133 = M2233 = 1 and hence that M = T. 
For each 3x3 matrix K and constant c let 



W c := Tr(K T E T (aI + cT)(EK)> = Tr(K' J 'AK) + cTr(K T E T T(EK)), 



(2.21) 



where I is the 4th order identity tensor. Throughout this paper we identify 3x3 matrices with 
9 dimensional vectors and 4th order tensors with 9x9 matrices following the usual convention: 



(1.1) ^1, (2,1)^2, (3,1)^3, (1,2)- 

(3.2) ^6, (1,3)^7, (2,3)^8, (3,3)- 

If we let A be the 9x9 matrix (4th order tensor) defined by 



A 0' 
A 
A 



4, (2,2)^5, 
9. 



(2.22) 



then W c can be written as 



hi 
hi 
hi 

h2 
h2 

fcl3 

h3 
.h3. 



(A + cM) 



hi 
hi 
hi 

kl2 

h2 
h2 

kl3 

h3 
.h3. 



(2.23) 



On the other hand, W c can be written as 

W c = Tr(K T E T L c (EK)}, 



(2.24) 
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where L c (x) is the 4th order tensor uI + cT and the constant c is chosen so that L c (x) is positive 
definite for all x. Since in the vector basis T takes the form 



T 



we have (in the vector basis) 
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(2.25) 



(2.26) 



Since Tr(K T E T T(EK)) is determined by the boundary data, the standard variational prin- 
ciple yields 



min (Tr(K W er(EK))) + cTr(K J W T(EK)) 
E = -VV — 

V = V° on dn 

min (Tr(K T E T L c (EK))). 
E = -VV _ 

V = V° on 



(2.27) 



We may rewrite this variational principle as 



mm 
E = -VV 

V = K T V° on dn 



(Tr(E T L c E)). 



(2.28) 



The constraints that (E) = I and V = K V u on 80, imply (E) = (EK) = K, and so we have 



W c > min (Tr(E T L c E)), 
(E)=K 

by relaxing the constraints. 

One can easily see that the minimum of the right hand side of (j2.29j) occurs when 



(2.29) 



and the minimum is 



L C E = n (a constant matrix) 
(Tr(E r L c E)) = Tr(K T (Lj 1 ) -1 



K). 



(2.30) 
(2.31) 
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Note that 

Tr(K T (L c : 1 ) _1 K) = hjUKY^ijklkkl- 
i,j,k,l 

So in the vector (rather than matrix) basis we get from (|2.23|) and (|2.2 

A + cM > (L" 1 ) -1 . 
By permuting rows and columns L c can be transformed to 



L,. 



a c c 
c a c 
c c a 



a —c 
—c a 



a —c 
—c a 



a —c 
— c a _ 

According to (I2.20P , A + cM takes (in that basis) the form 

Ai cMn22 cMn33 

cMn22 A 2 CM 2 233 * * 

cMn33 CM 22 33 A 3 

* A2 — cMn22 * 

-CM1122 Ai 

* * A3 -CM1133 

-cMn33 Ai 



A3 



-cAf 22 33 



-CM 22 33 A 2 



It is worth emphasizing that L c 1 takes the same block form as L c . By comparing 3 x 
on the left upper corner, f|2.33j) yields 



Ai 


CM1122 


cMn33~ 




O 


c 


c 


CM1122 


A 2 


CM 22 33 




C 


a 


c 


cMn33 


CM 22 33 


A3 




c 


c 


a 



-1 -1 



Note that the inverse of 



satisfy 



Thus we have 



ace 
c a c 
c c a 



is of the form 



abb 
b a b 
b b a 



where 



Kg + c) 

a = and 



a + 26 



(c-ct)(2c + ct)' 



1 



C7 + 2C 
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-1 -1 


"a /3 /3~ 


c 




c 




/3 a p 
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(7 




13 13 a 
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where 



a 



PM + (b)) 
(b) 



and j3 



(b) 



((b) - (a))(2(b) + (a)) 



(2.40) 



One can then see that as c is increasing and approaching a 2 , 
terms of eigenvalues) and approaching 



a c c 
c a c 
ecu 



-1 -1 



is decreasing (in 



d d d 
d d d 
d d d 



where 



3/: 



This shows in particular that 



a c c 
c a c 
c c a 



0"2 2o"2 + o\ 
1 -1 
) -02 



> 02- 



(2.41) 



1 1 1 
1 1 1 
1 1 1 



is positive-definite if c < o 2 re- 



gardless of the volume fraction f\. It then follows from (|2.36[) that 

-1 -1 
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0"2 



1 1 1 
1 1 1 
1 1 1 



(2.42) 



or 



Al CM1122 cMn33 
CM1122 A 2 CM 22 33 
cMn33 CM 22 33 h 
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a 


c 


c 


- cr 2 


1 


1 


1 


V < V- 


(( 


c 


a 


c 




1 


1 


1 






c 


c 


a 



-1 -1 



<72 



"1 


1 


f 
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Note that 





a 


c 


c 


« 


c 


a 


c 




c 


c 


a 



n -1 - 1 
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1 1 1" 




-1 


'a 13 13' 




1 1 1" 




7 5 5 


- 02 


1 1 1 
1 1 1 




l-( 


(3 a (3 
P a 


- o 2 


1 1 1 
1 1 1 


r- 


5 7 5 
5 5 7 



where 



5 (a - (72 + (3 - cr 2 ) 
P-02 



(3-a 2 



(2.43) 



(2.44) 



(2.45) 



(P - a)(2p - 2tr 2 + * - <r 2 ) 
Note that 5 blows up as c approaches a 2 , in phase 2. Thus to get a non-trivial bound we take 
v = [1 1 1] T and we have 



T 




7 (5 5 
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5 7 <5 
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5 5 7 
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3(7 + 25) 
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(<r+2c) 



3(72 — 02 ' 



(2.46) 



as c — > 02 where d is given by (|2.4ip . 

The left hand side of f|2.43j) approaches T as c —> 02 where 
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(2.47) 
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Thus we have the following upper bound: 



h < 



o\ + 2a 2 1 

(Tl - £J 2 1 + Cr 2 T' 



(2.48) 



This is generally not the tightest bound that can be obtained from (|2.33p as we have ignored 
the off-diagonal elements marked * in (|2.35p . A tighter bound (but one which is not so easily 
solved as an inequality on /1) is provided by the appropriate interval of volume fractions f\ 
given by the inequality (implied by (I2,33p ) 



det[A + a 2 M- (L" 1 )- 1 ] > 0, 



(2.49) 



in which, because L CT2 is singular, (L ) is defined as lim c ^. - 2 (L c } .In fact, one can show 
by tedious computations that in the basis in which L c is given by (I2.34p 



a a a 
a a a 
a a a 



-b b 



b -b 
-b b 



where 



h 



+ 



3/i 



b -b 
-b b 



11 
fi 



(2.50) 



(2.51) 



01 2cr 2 + G\ 

Note also that the translation tensor T is not the most general one: we could have taken T 
to be any symmetric fourth-order tensor satisfying the symmetries = —Tkju- We have not 
explored this freedom in the choice T, which could lead to even tighter bounds. 



2.2 Special boundary conditions 

If the affine Dirichlet boundary conditions are prescribed, namely, V° = x on <9f2, then (E) = I, 
and hence 

A = (E T aB) = (E T ){aB) = Ia D (E) = a D , (2.52) 
where ajj is the Dirichlet tensor. Since M1122 = -M1133 = -W2233 = 1 in this case, we have 

T = Tr[(a D -a 2 I)- 1 ], 



and hence 



fi< 



0\ + 2cr 2 



CJi - <72 1 + cr 2 Tr[(<X£> - cr 2 I) x \ 

This upper bound is exactly the same as the one obtained in |18j - 
The inequality f|2.33j) yields not only (|2.36p but also 



(2.53) 
(2.54) 



A 2 — cMn22 

-CM1122 Ai 

A3 -cMn33 

-CM1133 Ai 
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—c 






—c 
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—c 
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a 



(2.55) 
(2.56) 



S 



and 



A 3 -cM 22 33 
-cM 22 33 A 2 





a —c 






—c a 



(2.57) 



When V° = x on <9f2 these inequalities do not yield a better bound than (|2.48|) . To see this, 
observe that with these boundary conditions (|2.55p becomes 



A 2 — c 




a —c 






— c Ai 




—c a 



-1 



and hence 



A 2 -c 
— c Ai 



0"2 



1 -1 
-1 1 



> 



a — c 
— c cr 



-i 



-i 



1 -1 
-1 1 



We then obtain 



1 1 1 
+ < 



Ai — cr 2 ' A 2 — <7 2 ! 



where 



^ 2<T 2 <T 2 + (Tl / 

So far we have shown that (I2.55f) . (|2.56j) and (I2.57P yield 



/2 /l 



1 



+ 



1 



+ 



1 



< 



Ai — cj 2 A 2 — cr 2 A3 — cr 2 7— <r 2 



which in turn yields 



Note that 



/i< 



(Tl + <7 2 



1 



0-1 - ct 2 1 + fcr 2 Tr[(<r £) - a 2 I)- 1 ] 

3 



Tt[(o- D - ^I)- 1 ] > 



(Tl - (T 2 



Using this fact, one can easily show that 

(Tl + 2(T 2 



< 



(Tl + <r 2 



1 + a 2 Tr[(a D - a^I)- 1 ] ~ 1 + |a 2 Tr[(cr D - a 2 I)-l] 



(2.58) 



(2.59) 

(2.60) 
(2.61) 

(2.62) 

(2.63) 
(2.64) 

(2.65) 



Thus the bound in (|2.54j) is better than the one in (12.63p . They are the same bounds only when 

h = 1. 



3 Attainability of the determinant bound (12.491) 

Suppose that the bound in (|2,29p is attained for a field E. Then one can show (see for 
more details in the analogous two-dimensional case) that 

L CT2 E = (L; 2 1 )- 1 (E), (3.1) 

where (L~ ) , given by (1230]) . is defined as lim c _ ) . cr2 (Lj 1 ) -1 . We now show that the attain- 
ability condition (13.11) holds if and only if 

L ff2 E = U, (3.2) 
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for some constant matrix U. The only if part is trivial. The converse is clear heuristically: 
Since E = L^U and U is constant, (E) = (L~£)U and hence U = (L- 1 ) -1 (E), therefore (gj} 
holds. To show the converse rigorously, we use the basis in which L c is given by (|2.34j) . In that 
basis, 

E = [-En E22 E33 E21 E\2 E31 E13 E32 E23] . (3.3) 

Let an d denote L CT2 in phase 1 and 2, respectively, and likewise E*- 1 ) and E 1 - 2 ). Since 

(2) 

U is in the image of L^ 2 , U takes the form 

U=[aaap 1 -ft ft -ft ft -ft] T . (3.4) 

Thus E^ 2 ) takes the form 

E( 2 )=E + E ± , (3.5) 

where 



Er 



a a a 01 01 02 02 03 03 

3(72 3ff2 3(72 2(72 2(72 2(72 2(72 2(7 2 2<r 2 

(2) 



T 



(3.6) 



and Ej_ is in the kernel of L^2 , i-e., Ej_ is symmetric and Tr Ej_ = 0. Since Jph ase 2 ^j. * s a ^ so 
symmetric and of trace 0, we can see from the form (|2.50p of (L^ 1 )^ 1 that 



j)- 1 [ E ± = 0. (3.7) 
Jnhase 2 



(L 

(72 

'phase 2 
Thus we obtain 

(L- 1 )- 1 ^) = (L- 1 )- 1 (/i(L«)- 1 U + / 2 E ) = U, (3.8) 

where the last equality can be shown by an elementary but tedious computation. 

The condition (|3.2p implies in particular that E is constant in phase 1. We now show that 
this condition alone guarantees that the upper bound is attained. 

Suppose that E is constant in phase 1. After coordinate changes and taking linear combi- 
nations of potentials if necessary, we may assume that 

E = VV = 1 in phase I. (3.9) 

We claim that if (|3.9h holds, then there exists a potential ip in phase 2 such that 

V = Vip and Aip = ° X + 2 °" 2 in phase 2. (3.10) 

To show (|3,10p let D be a connected component of phase 1 and J7o an open subset of SI 
containing D such that 0,q \ D is phase 2 and connected. It is known (see section 2.5 of [4j) 
that for j = 1, 2, 3, there is a harmonic function Hj and a potential <pj such that 

Vj(x) = iZ^x) + / ^' (y) da(y), xel] . (3.11) 

JdD 47r l x - y| 



Furthermore, we have 



If (EH) holds, then 



_ o\ — 02 dVj 
Vj ~ 02 dn 



on dD. (3.12) 



— Zl — —rij on dD, (3.13) 

0"2 
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and hence 



Thus we have 



cri - a 2 d f 1 
<?2 dxj J D 4vr|x - y| 



Vj (x) = Hj (x) - — !— — - —/ ,_,„ .,, dy, xeH . (3.14) 



flj-(x) = ^ + e,- + A / — -i -dy, xefl, (3.15) 

cr 2 y D 4vr|x - y| 

for a constant Cj. Let H = (Hi, H2, H3) and c = (01,02,03). Then, 

H(x) = v(^|x| 2 +c-x+^^ / — -i -dyV x£D. (3.16) 

V 2 o- 2 y jD 4vr|x-y| J 

In particular, we have V x H = in D. Further we have V H = ?±±2£± - m d. Since all 
components of V x H and V • H are harmonic in Qq, we have V x H = and V • H 



a ■> 



VLq. It then follows from (|3.14p that V x V = and V • V = ^ in Qq \ D, and hence in phase 
2 since the phase 2 is connected. 

This does not quite guarantee existence of a potential tp such that V = VV> since phase 2 is 
not simply connected. But, V can be written as V = H + S where H is harmonic in Q and S 
is harmonic in R 3 \ {phase 1} and S(x) — > as jx| — > 00. What we have shown above is that 
V x H = in Q, and V x S = in I 3 \ {phase 1}. It implies that there are potentials and 
ip2 such that H = V^i i n ^ an d S = in R 3 \ {phase 1}. Thus we have (|3,10p . 

The relation (|3.10p implies that E is symmetric and TrE is constant in phase 2. Therefore 
L CT2 E is constant in phase 2. Because of the continuity of uE along the interface, it follows that 
L CT2 E is constant in £1 and hence the determinant bound ()2.49p is attained. 



4 Lower bounds 



In this section we derive lower bounds for f\. We still use the translation method. But as 
there is no null-Lagrangian for current fields (see below) we use a quasiconvex quadratic form, 
namely the extremal one discussed on page 546 of [17] . 

Let Ah, A s and A a be orthogonal projections onto the space of 3 x 3 matrices proportional 
to I, symmetric and trace free, and antisymmetric, respectively, namely 

A ft P = ^Tr(P)I, 

A S P = i(P + P T )-A h P, 

A a P = ^(P-P T ), 

for any 3x3 matrix P. Define a 4th order tensor T' by 

T' := 2A S - A h , so that T'P = P + P T - Tr(P)I. (4.1) 

This tensor is positive semidefinite on the set of 3 x 3 rank two matrix valued fields, and isotropic 
in the sense that 

T'(R r PR) = R T (T'P)R, (4.2) 

for any unitary transform R. Due to this isotropy it suffices, to check for positive semidefinite- 
ness on rank two matrices, to consider matrices of the form 



Pn P12 P13 
P21 P22 P23 




(4.3) 
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and in that case 

Tr(P T T'P) = ( Pll - V22? + (P12 + P2i? + P 2 13 + vis (4-4) 

is clearly non-negative and zero only when p\\ = P22, P12 = — P21, and P13 = P23 = 0. More 
generally, Tr(P T T'P) will be zero if and only if the rank two matrix P is of the form 

Pij = a (kikj - Sijk ■ k) + f3 e ijm k m , (4.5) 

for some constants «o an d fio, and for some vector k which is a vector such that Pk = 0. 
Here eum is the completely antisymmetric Levi-Civita symbol, taking value 1 if klm is an even 
permutation of 123, —1 if it is an odd permutation, and zero otherwise. 

If J(x) is periodic, with unit cell T, and has zero divergence then its Fourier components 
J(k) satisfy k • J(k) = and hence for k 7^ are matrices of rank two or less. As observed by 
Murat and Tartar |22[ [20"1 [24"] , it then follows by a straightforward application of Plancherel's 
theorem that the quadratic form associated with T' is quasiconvex on divergence free periodic 
fields: 

(Tr(J T T'J)) T > Tr((J T ) T T / (J) T ), (4.6) 

(where the angular brackets (-)x denote volume averages over T) with equality if and only if 
the Fourier components of J have the form (14. 5p . i.e. 



to 9 2 a _, . df3 , A _ S 

lkl = Jkl + 5^ " 6kAa + eklm d^ (47) 

for some scalar potentials a(x) and /3(x), where the are the elements of a constant matrix 
J°. Note that J given by (|4.7p is divergence free, has J° as its average over T, and 

T'J = T'J° + 2VVa (4.8) 

is the gradient of a potential. 

Incidentally, there is no symmetric fourth order tensor T', with elements T^ M , such that T' 
is a null Lagrangian for currents. If there were the quadratic form Tr(P T T'P) would be zero 
for all 3 x 3 rank two matrix valued fields. In particular, being zero for all matrices of the form 
(|4.3p implies T[- kl = unless i = 3 or k = 3. Similarly, by considering matrices where the first 
or second row of P is zero we deduce that T[- kl = unless i = 1 or k = 1, and T[- ki = unless 
i = 2 or k = 2. It then follows that all elements of T' must be zero. 

Let ji, J2 and J3 be the currents corresponding to the measurements ji ■ n = q\, J2 ■ n = #2 
and j3 • n = on <9Q. Let J = [ji J2 J3] and suppose that 

(J) = I. (4.9) 

Define the response matrix A' by 

A' := (3 T a- 1 J), (4.10) 
which is computable from the boundary data: 

A '=wI jT( - vv) =wL- (nJ)Tv - (4 ' n) 

For a positive constant c, define 

h' c := a' 1 ! - cT' = (a' 1 + c)A h + (cr^ 1 - 2c)A s + (T~ l k a . (4.12) 

We assume that 

c < (4.13) 
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so that \J C is positive definite. It is helpful to recall that a 2 < a±. 
Define 



where q = [q-y q 2 qs}. Let 



Then we have 



W':= min Tr(J J L' r J), (4.14) 

v-j=o, (J)=i, 

nJ_— — q on 



g : = min Tr(J T T'J). (4.15) 

3 v-j=o, (j>=i, \— — ' ^ ' 

nJ_— — q on 



W'< min Tr(J T <T _1 J) - min cTr(J T T'J) = Tr A' - eg < Tr A' - cg~ , (4.16) 
c _ v-j=o, (j>=i, s - -' v-j=o, (J>=i, \— — / 3 - » ) v y 

nJ_— — q on OH nJ_— — q on dfl 

where g~ is a lower bound on 5. Lower bounds on g are easily obtained. Choose a cube T 
which contains £1 and construct a divergence free field J in T \ f2 such that J satisfies periodic 
boundary conditions on <9Y and nJ = — q on dil. Then from (I4.6P the inequality g~ < g holds 
with 

g~ = - {Tr[(pl+ (1 -p)(J)t\q) T T'(pI + (1 -p)(J) r \ n )} ~ U - p) Tr[(J T T'j) T ^]} , (4-17) 

where p is the volume fraction the body £1 occupies within the cube T (and the angular brackets 
(■)x\n denote volume averages over T\0). This bound is not necessarily useful unless J is close 
to being of the form (|4.7p . There is a wide variety of fluxes q for which g is exactly computable. 
Let smooth potentials a and j3 be chosen in the neighborhood of the boundary of f2 and set 



« = + 5 M Aa + e k i m - — . (4.18) 



8 2 a . . dp 
kl ^d^ l - 6klAa + eklm d^r 

Then we can extend a and f3 to all of T with periodic boundary conditions on dT, and thus 
obtain a field J given by (|4.7p satisfying nJ = — q on d£l. This field attains equality in (|4.6p : 

P Tr[(J T T'J)] + (1 -p)Tr[(J T T'J) r \ n = Tr[(J°) T T'J ]. (4.19) 

The first term in this equation cannot be reduced, while maintaining the inequality ()4.6p . by 
varying J inside f2 while keeping nJ = — q on d£l (which ensures (J) remains equal to I). Using 
(14.8P and integration by parts this establishes that 

g = TV[(J T T'J)] = -L / -[x T T'J° + 2Va] • q. (4.20) 
\M\ Jan 

In particular, if we take special Neumann conditions q = — n on <9f2, then (|4.18p is satisfied 
with J° = I and a = j3 = 0, and the first equality in (|4.20p implies g = —3. It is an open 
question as to whether given any fluxes q\, q 2 and q% (with no net flux through dQ) one can 
find a, j3 and J° such that (|4.18p holds, however the following calculation suggests it might be 
possible. Suppose that there is a portion of dQ which is flat with normal n = (0, 0, 1) T . Then 
on this flat section (|4.18p implies 

o da 3 8[3 _ da 3 80 Q d 2 a 8 2 a 

qi -- J ^-^[-^ q2 -- J ^~d^ + ^ q3 ~- J33+ d4 + ^f' (421) 

where 03 = 803/8x3. Thus possible values of a, 03, and ft on this flat interface may be found 
by solving Poisson's equations 

<9^a + c^a = +j0 d 2 a 3 | 8 2 a 3 _ ggj 8q 2 cPp (?P_ _ dq^ _ 8qi_ ^ ^ 
8x\ 8x\ qs 33 ' dx\ 8x\ 8x\ dx 2 8x\ 8x\ 8x\ 8x 2 
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Now we have from (|4.12p 



W' c > min Tr^er" 1 + c)A h + (a" 1 - 2c) A s + a^AJJ) 
> mmTr^iia' 1 + c)A h ]J), 



(4.23) 



by relaxing constraints. Observe that 

(J^Ka- 1 + c)A fc ]J> = ^((a- 1 +c)(TVJ) 2 } = ^J_+£ 
It thus follows from Jensen's inequality that 



(TrJr + 



2 , °2 + c 



phasel 



3\n\ 



(TrJ) 2 



phase2 



CTi 1 + C / 1 



101 



3/i 



2 -1 , /i 

(T 2 + C / 1 



phasel / •J/2 



Tr J 



phase2 



We also have from the constraint (J) = I that 



1 



— / i v. i • -!- 



\n\ 



phasel 



n\ 



TrJ = 3. 



phase2 



By minimizing the righthand side of (|4.24j) under the constraint f|4.25|) we get 

3 



W' > min TV^ffcr" 1 + c)AJJ) 
(j)=i 



A 



(4.24) 



(4.25) 



(4.26) 



for all c < (J 1 /2. Letting c — >■ cj 1 /2, we have from (|4.16|) that 



3<x, 



, 2/ 2 
r T TT^T — 



which yields 



fx = i - h > i 



2<7l + <7 2 



2(T 9 
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2cri Tr A' - g" 



(4.27) 



(4.28) 



2(<Tl - <7 2 ) 

With the special Neumann conditions q = — n on d£l, we have (J) = I and 

A' = {J T a~ l J) = (J T )(a- l J) = Ia N l (J) = * N \ (4.29) 
where is the Neumann tensor. Since g = —3 in this case, we can take g~ = —3 giving 



fx>l 



2<7i + cr 2 



1 



9 



2a! Tr(<7^) + 3 



2(<7i - CT 2 

As for the upper bound, a lower bound was obtained in |18| : 

2fTi + cr 2 1 



(4.30) 



fx>l 



(Tl - cr 2 <7i Tr[(<7il - o-jv) x ] - 1 ' 



(4.31) 



The bound in (|4.30[) coincides with that in (|4.3ip when = <7jvI. We could have recovered 
the bound (|4.3ip , if we had defined 

W' := min Tr(K T J T U(JK)}, g := min Tr(K T J T T'(JK)), (4.32) 

c v-j=o, (J>=i, x — cy — " VJ=0, (J>=I, \ — V— //' v > 

nJ — — q on dQ, nJ — — q on dfl 
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and then optimized over the choice of K, in a similar manner as was done for the upper bound. 
This optimization is straightforward with the Neumann conditions q = — n on dCl, but in 
general might be difficult as it is not clear if the dependence of g on K is quadratic. A way 
around this difficulty is to fix (independent of K) a J in T \ SI such that J satisfies periodic 
boundary conditions on <9T and nJ = — q on <9$7 and to take 

g~ = X - {Tr[(pK + (1 - p)(J) T ^K) T T'( P K + (1 - p)(J) T \ n K)] - (1 - p) Tr[((JK) T T'(JK)) T \ n ]} , 

(4.33) 

which then depends quadratically on K. However in the end one is still faced with the problem 
of choosing J in T \ Q, to optimize the bound. The treatment given here can be generalized in 
other ways too, notably by allowing for translations T' which are not isotropic. 
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